Theory of elastic anomalies at the 7—7-0! phase transition in solid Ce 
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Starting from a model of 4f-electron generated quadrupolar densities on a compressible fee lattice, 
the elastic anomalies at the Fm3m — > Pa3 phase transition are studied by means of analytical theory. 
The model is taken as representative for the 7 — a phase transition in Ce. The coupling of the (linear) 
lattice displacements to the square of the quadrupolar orientational density fluctuations renormalizes 
the elastic constants. The condensation of the quadrupolar densities into the orientationally ordered 
Pa3 structure is studied as function of temperature and pressure. Precursor effects of the transition 
lead to an anomalous softening of the elastic constant en while C44 exhibits no such softening. The 
theoretical results are in excellent qualitative agreement with pressure experiments on the elastic 
constants (equivalently on sound velocities) at the 7 — a transition in Ce. Lattice dynamical analogies 
in theory and striking similarities in experimental results with the Fm3m — > Pa3 transition in Ceo 
fullerite are discussed. 

PACS numbers: 61.50.Ks, 62.20.de, 64.70.K-, 63.20.K- 
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I. INTRODUCTION 



Cerium, the first of the lanthanide series (rare earths) 
with an inner 4f-electron exhibits a broad range of un- 
usual electronic, magnetic and structural properties in 
the solid state [l|. Most intriguing is the transforma- 
tion under pressure Q between two face-centered-cubic 
(fee) phases called 7-Ce and a-Ce. While 7-Ce is sta- 
ble at room temperature (T) and ambient pressure (P), 
it transforms to a-Ce at F w 0.8 GPa. This appar- 
ent isostructural phase transition 7 — > a is accompa- 
nied by a ~15% volume reduction and a change in mag- 
netism. While in the 7-phase the magnetic susceptibil- 
ity x(T) follows a Curie Weiss law which is attributed 
to the magnetic moment of the 4f-electron, it is a quasi 
T-independent in the a-phase and reminiscent of Pauli 
paramagnetism Q. 

The isostructural nature of the 7 — > a transition is 
at variance with the Landau theory of phase transitions 
[3] which would require a change of space group symme- 
try. Since the establishing of the fee structure for both 
phases Q various microscopic theories have been pro- 
posed. Nowadays two theoretical models are well repre- 
sented in the literature. 

A Mott-transition scenario for 4f-states is based on the 
concept of a localized non bonding state of the 4f-electron 
in 7-Ce and an extended metallic bonding state in the a- 
phase @. The increased binding capability due to the 
metallic f state is considered then as the reason for the 
marked volume reduction in the a-phase @. In the fol- 
lowing this model has been supported by band structure 
calculations |7[ . On the other hand more recent inelastic 
neutron scattering experiments on a-Ce found that the 
magnetic form factor is quite different [§J from the one 
calculated [1] under the assumption that the 4f states are 
itinerant. The absence of extended 4f states in the solid 
is also supported by many-electron calculations on the 



Ce dimer [1CH12| where it is found that the 4f-electrons 
do not participate in the chemical bond. 

The second widely recognized theoretical model is 
based on spin fluctuations. These are the driving mecha- 
nism in the Kondo-volume-collapse (KVC) model of the 
7 — > a transition [l3|, E3] ■ In first approximation the 4f- 
electron stays localized in both phases, however in the 
a-phase the hybridization between the 4f-electrons and 
the conduction electrons (Anderson impurity Hamilto- 
nian, see e.g. Ref. Il5r ) is much more intense than in the 
7-phase. This interaction leads to a screening of the 4f 
local moments in the energy ground state. The volume 
dependence of the hybridization coupling or equivalently 
of the Kondo-temperature Tk is used to obtain the equa- 
tion of state. The interpretation of the 7 — > a transition 
in Ce within the KVC scenario [HI has found stron g su p- 
port by the analysis of electron spectroscopy data [16| ■ 

Notwithstanding much effort, both the Mott transition 
and the KVC model remain under debate [17| and there 
is no consensus on the driving mechanism of the phase 
transition fisl . . A common feature is the absence of 
any symmetry breaking at the transition as would be 
required by Landau theory [4j. On the basis of ther- 
modynamical data it has been suggested [2(j that the 
7 — > a transformation is in fact a first order phase tran- 
sition which becomes of second order beyond a tricritical 
point. Therefore, a-Ce should have lower symmetry than 
7-Ce. However a distorted lattice has been discarded by 
x-ray diffraction experiments pH l22l|. 

A mechanism of symmetry lowering without lattice dis- 
tortion (positions of the Ce nuclei fee in the 7 and a 
phase) has been suggested by the present authors [23l — 
|25| . The main concepts are the support of quadrupolar 
electronic charge-density fluctuations by the 4f-electrons 
in the 7-phase on a compressible fee lattice and the col- 
lective orientational ordering of the quadrupolar densities 
in the a phase (23W25I ] on four simple cubic (sc) sublat- 
tices. The space group symmetry lowering at the transi- 



2 



tion 7 — > a is Fm3m —> Pa3, the local density symmetry 
is Sq. It is accompanied by a uniform lattice contraction 
so that the fee structure is conserved. Although this lat- 
ter aspect is "isostructural" , the symmetry lowering of 
the electronic structure is fully consistent with Landau 
theory. In reciprocal space the active electronic density 
mode condenses at the X point of the Brillouin zone. 

A phase transition Fm3m — » Pa3 occurs in solid Cgo 
(fullerite) [26, 27]. Here the symmetry lowering is due to 
the orientational ordering of the Ceo molecules. At the 
transition the cubic lattice constant contracts discontin- 
uously [H, , the center of mass points of the molecules 
still occupy an fee lattice. From the latter point of view, 
the phase transition would be "isostructural" . However, 
the molecular order on four sc sublattices entails an order 
of the constituent C atoms and hence the Pa3 structure 
leads to characteristic reflections in x-ray [26| and neu- 
tron scattering [27j experiments. 

The problem of measuring the Pa3 structure directly 
in Ce is complicated by the fact that the Ce nuclei re- 
main on an fee lattice. Although the quadrupolar charge 
density ordering could be detectable by synchrotron ra- 
diation, the intensity of the additional Pa3 reflections is 
likely to be too weak and has not yet been uncovered 
30]. On the other hand there have been two recent ex- 
perimental results which are very specific and point to 
the relevance of the quadrupolar ordering at the 7 — > a 
transition. 

Time-differential perturbed angular correlation (TD- 
PAC) experiments in solid Ce in a pressure range up to 8 
GPa detect an appreciable electric field gradient (EFG) 
in a-Ce which is almost four times larger than in the cu- 
bic 7 phase and close to values in the noncubic phases 
a' and a" [31]. This finding rules out the Fm3m sym- 
metry in a-Ce and evidences in support of the antifcrro- 
quadrupolar order suggested by theory [23l425l |. 

The other experimental support of the antiferro- 
quadrupolar ordering is more indirect but very specific 
too. Phonon dispersion measurements by inelastic x- 
ray scattering on elemental Ce across the 7 — > a transi- 
tion reveal strong changes in the dispersion shape [32] ]. 
In particular a pronounced softening of certain phonon 
branches is found in the a-phase toward the X point of 
the Brillouin zone. 

Given the above mentioned fact that the discussion on 
the validity of the Mott transition or KVC scenario re- 
mains open, there have been in the last decade increased 
experimental efforts in investi gating the lattice dynamics 
at the phase transition 7 — > a [2lL I22I I32I HH • From these 
studies it results that the lattice vibrations play an im- 
portant role at the 7 — > a phase transition and that the 
inclusion of lattice dynamics in the theoretical descrip- 
tion is of paramount importance. In this respect we recall 
that anomalous elastic behavior at the 7 — > a transition 
was originally discovered in ultrasound experiments [3~ij . 
From inelastic neutron scattering experiments [35[ per- 
formed on single crystal of 7-Ce it was concluded that 
premonitory effect of the 7 — a transition are present 



in the phonon dispersion curves of 7-Ce at room tem- 
perature. However, until recently the relevance of these 
results has not been sufficiently appreciated, while all re- 
search efforts were concentrated on the electronic prop- 
erties. 

The anomalies of the elastic properties at the 7 ^ a 
transition in Ce [34,21,33] have their counterpart at the 
fee— >sc transition in Cgo fullerite [36l |37| . In particular 
there is a striking resemblance of the pressure dependence 
of the bulk modulus in Ce with corresponding experimen- 
tal results in C 60 [13 ■ 

The content of the paper is as follows. In Sec. |TT] we 
recall the main features of the model of interacting 4f- 
electron quadrupolar densities on a compressible fee lat- 
tice. The Hamiltonian and the resulting condensation 
scheme from the disordered 7-phase to the quadrupolar 
ordered a-phase are described. The coupling to the lat- 
tice is included. Next (Sec. IIII| we extend the model by 
adding an external hydrostatic pressure. Various quan- 
tities that characterize the phase transition (order pa- 
rameter, susceptibility and correlation length) are de- 
rived as functions of temperature and pressure. In Sec. 
IIVI we calculate the dynamic displacement-displacement 
correlation function, taking into account the coupling of 
the quadrupolar electron densities to the crystal lattice. 
In the static limit the inversion of the displacement- 
displacement correlation function tensor leads to the elas- 
tic constants Cn and C44. While in cn the coupling 
to quadrupolar electron density fluctuations leads to re- 
markable anomalies at the 7 — a phase transition, the 
coupling is absent in C44. The results of the theory are 
compared with experiments. Conclusions of the paper 
are presented in Sec. [V] 

II. MODEL 

Here we will remind the model Hamiltonian which de- 
scribes electronic quadrupolar charge density fluctuations 
on a compressible fee lattice. In previous work we have 
first treated quadrupolar charge fluctuations due to solely 
4f-electrons [23|. Later on we have taken into account 
conduction electrons 1241 as well as f- and d-electron in- 
tra site correlations [25| . Since these extensions can be 
cast into an effective Hamiltonian with essentially the 
same structure as the one originally studied [Hj], we will 
restrict ourselves here to 4f-electrons as a generic case. 

We consider N Ce atoms located on a non rigid fee lat- 
tice with nuclear positions R(n) — X(n) + u(n), n = I, 
2, ... , N. Here {X(n)} denote the equilibrium positions 
and {u(n)} the displacements due to lattice vibrations. 
We assume that core and valence electrons follow adi- 
abatically the nuclear displacements. The Hamiltonian 
comprises three parts: 

H = + P (L) + H {e - L \ (1) 

where stands for the 4f-electrons, H ' L > for the lat- 
tice and iP e ' L ) for the coupling of both. 
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Explicitly one has 

=Uo + U QQ , 



(2) 



where Uq is the single particle potential and Uqq 
the quadrupole-quadrupole interaction on the rigid lat- 
tice. We assume that the 4f-electrons have coordinates 
{X(n) + f(n)} and are localized on spheres centered at 
jX(ri)} and with radii {f(n)}, \r(ft)\ — rf =1.378 a.u. 
[24j [. The single particle potential is due to the cubic 
crystal field in presence of spin-orbit coupling. One ob- 
tains 



U 



EE 



(3) 



where and \i)f L are the eigenvalues and eigenstates, 
i = 1, 2, ... , 14. The angular part of the wave function 
of the 4f-electron (angular momentum quantum number 
I = 3) at site n is described by the 2(21 + 1) spin orbitals 
(i\n), n = @(n), where 8 = (&,(/>)■ For details see Refs. 
I23l425l . The electronic quadrupolar interaction is given 

by 



U, 



QQ 



p^(n)v a? (n-H')p^(ft), (4) 



up 



where Pq {ft), oi(fi) =1, 2, 3 are quadrupolar density 
operators and v a ^(ft — ft') interaction matrix elements 
between nearest neighbors on the fee lattice. The density 
operators are defined by 



PqH™) = E C Q w')N)nvi|n! 



with 



(ij)= / dn{i\n)S^(n){n\j). 



(5) 



(6) 



Here we restrict ourselves to the three quadrupolar func- 
tions S'W(n), a = 1, 2, 3, which transform as the irre- 
ducible representation Ti g of the cubic point group m3m 
(Oh)- These symmetry adapted functions (SAF's) are 
linear combinations of spherical harmonics belonging to 
the manifold I = 2 and are tabulated in Ref. [3|| They 
transform as the Cartesian components yz, zx, xy for 
a = 1, 2, 3, respectively. Defining Fourier transforms 



1 



E< 



-iq-X(n) (a) 



Pq in), 



(7) 



v a0 (q) = e- lf( * (S,) -* (S) Wn - ft), (8) 



we obtain 



U QQ = T^EE^®^®^ - ^' 



(9) 




FIG. 1: Triple-<f antiferroquadrupolar (Pa3) structure of a- 
Ce proposed earlier [53]. Quadrupoles represent the I — 2 
valence electron (4/ + 5d6s 2 ) charge density distribution. 



with [25 J 



v(q) 



'■yC yz + a(C zx +C xu ) -pS^y 

-PS ay jC zx + a(C XH +C y z ) 

-ps zx -PS yz 7 C X 



-ps zm 
-ps yz 

■ + a(C yz +C z!c ) 



(10) 



Here 7, /3 and a are quadrupole-quadrupole interaction 
coefficients on the fee lattice. The quantities Cy = 
cos{qia/2)cos(qja/2), and SV,- = sin(gja/2) sm{qja/2), 
i = x,y,z, a is the cubic lattice constant, account for 
the fee lattice structure. At the X point of the Brillouin 
zone for T£, the interaction matrix v(q) becomes diago- 
nal. The star of q x has the arms q£ = (2ir/a)(l, 0, 0), 
qf = (27r/o)(0,l,0) ) gf = (27r/o)(0, 0, 1). For q = qf 

one has V22((lx) — v 33{Qx) — ~^J> an d similarly for 
q* and 5^ with permutation of indices (act). Hence 
the quadrupolar interaction becomes attractive with the 
largest twofold degenerate eigenvalue X x + = at each 

arm of q~f . 

As we have shown previously [23 25] this interaction 
is compatible with a symmetry lowering Fmim — > Pa3, 
characterized by an orientational order of the quadrupo- 
lar densities on four different sc sublattices (see Fig. [T|) . 
The corresponding condensation scheme reads: 



,(3)< 



(l)e 
= PQ 



(qX) = p% )e (q?)=pVN^0, (11a) 



(2)e/-.X\ _ „(3)e/-i?s _ „(!)« 



(<ff) = 0. (lib) 



9 a P 



Here the superscript e stands for a thermal expectation 
value and p is the order parameter amplitude. Since the 
three arms of q%? are involved one speaks of a triple-g an- 
tiferroquadrupolar order (AFQ). Notice that a conden- 
sation scheme similar to Eqs. (jllap . (lllb[) also holds for 
the phase transitions with the symmetry Fm3m — > Pa3 
in NaC>2 [39| and in Ceo-fullerite [4(|. Although icosa- 
hedral symmetry of the Cgo molecule implies that not 



4 



quadrupoles but higher order multipoles (I = 6, 10, 
...) determine the orientational interactions, the corre- 
sponding SAF's transform as irreducible representations 
of symmetry Ti g of the cubic point group. 

A compelling mathematical reason to consider the 
symmetry lowering Fm3m — > Pai as a candidate for 
the isostructural phase transition in Ce is the fact that 
it leads to a lattice contraction while the center of mass 
points still occupy an fee lattice. In addition the tran- 
sition is of first order. Indeed the number of symmetry 
elements is thereby reduced by a factor 3 (from 48 to 
16), which implies the existence of a third order cubic 
invariant in free energy jij. 

The lattice dynamics is described by the phonon 
Hamiltonian 



H {L) = K + U, 



(12) 



where K is the kinetic energy and U the potential energy 
in harmonic approximation. In Fourier space one has 



1 



(13a) 



P = ? EE M ii(<?>l (13b) 



= x, y, z. The displacements Ui{q) and the conjugate 
moments Pi(q) are related to the variables in real space 
by 



Ui(n) 



Pi(n) 




Ui(q)e 



iq-X(n) 



iq-X(n) 



(14a) 



(14b) 



displacements and for quadrupolar fluctuations near the 
X point of BZ we have 



^'^EEEE'^'^ 



P-Q)Pq (£) 



(17) 



Here the prime on the sum over a indicates the following 
restrictions: for i = x, a = 2, 3; i = y, a = 3, 1; i = z, 
a = 1,2. The coupling matrix is given by 



v'i, aa (q,p) = 



= J2(Q- aa {h) cos(p ■ X(h)). 

Nm i ~T 

h 



(18) 



Here X(h) = X(n') — X(n), ft 1 refers to the 12 nearest 
neighbors of n on the fee lattice, v\ aa {n' ~ n) is the first 
order derivative with respect to the lattice displacement 
component i of the quadrupole-quadrupole interaction 
v aa (ft' — ft) . The structure of the fee lattice implies that 

^,33(^1) = ^,22(^2) = v' itll (h 3 ) = A, where for i = x 
or y, X(hi) = a/2(l, 1,0); for j = z or x, X(fi2) = 
a/2(l, 0, 1); for i = z or y, X{h 3 ) = a/2(0, 1, 1). Within 
the quadrupolar model we obtain that A < 0. Carrying 
out the summation over h and exploiting the symmetry 
of the lattice, we rewrite Eq. (fTT|) as 



zA<2 ^ — r ^ — > ^ — "v ^ — V 

8 



qiUi{q) 



p 
» 



(-P-Q)Pq (f) 



(19) 



Under the proviso of the summation restriction over a 
we define for i — x, y or z 



" w (?) = E J2'pq ) (-p + ftp™ ® . 



(20) 



where m is the atomic mass. One has the usual commu- 
tation rules [it, u] = \jp,p] = and 



(15) 



In the long wavelength limit the dynamical matrix M(q) 
is given by 



M(q) 



9x c 



4m 

,.J, 



+ c 4 
1 r o 



'A, + d 

9z<7 B ( c 12 



+ ^4) 



4*4l +(Ql+Q 2 y )< 



and hence 




^^^(g)^)(-q). (21) 



Only longitudinal lattice displacements or equivalently 
longitudinal lattice strains en occur on the right hand 
side of Eq. (|19p . Since in addition qiUi for i — x.y and 
z occur on the same footing, the coupling H^ e ' L ^ leads 
to a striction of the lattice with conservation of cubic 
symmetry 23] . Notice that there is no coupling to trans- 
verse lattice displacement waves or equivalently to shear 
strains ey, i ^ j, in H^ e ' L \ In the next section we will 
investigate the influence of the quadrupolar density fluc- 
(lg)tuations on the low frequency lattice dynamics. 



where 



are the bare elastic constants in absence of 



coupling to the quadrupolar electronic fluctuations. 

The coupling between quadrupolar charge density fluc- 
tuations and lattice dynamics has been derived previ- 
ously 23j. In the long wavelength regime for the lattice 



III. PHASE TRANSITION 

In order to study the tripple-<f antiferroquadrupolar 
phase transition on a compressible lattice as a function of 
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TABLE I: Parameters x {2) , x (3) , 7, A, B, C, C , D calculated 
from Ref . I23I. Pi and kl estimated from experiment [3^ ]. 



x {2) = 23.34 x 10" 
47 = 3491 K 
B = 8813 K 
C = 103692 K 
a 7 = 5.169 A 

kl = (37.0 GPa)" 



x (3) = -1.30 x 10" 
A = 445 K/A 
C = 103964 K 
D = 216700 K 
a a = 4.857 A 
Pi = 0.8 GPa 



temperature and pressure, we first recall some concepts 
of the underlying free energy. The phase transition is 
of first order. We calculate the order parameter as a 
function of temperature and pressure. We show that the 
transition temperature increases linearly with pressure. 
Finally we study the order parameter susceptibility in 
the disordered and the ordered phase. 

Taking into account the quadrupole-quadrupole inter- 
action on a rigid fee lattice we have written the Helmholtz 
free energy as a Landau expansion in terms of the order 
parameter amplitude p [23]: 



Fqq/N = F /N + Ap 2 + Bp 3 + Cp A . 



(22) 



Here Fq is the free energy of the disordered phase which 
has to be calculated with the cubic crystal field Uq, 
Eq. ([3]). The coefficients of the order parameter terms 
are 



A = - 
2 

B = 



C 



T 

x {2) 4 T 



(23) 
(24) 



T 



9(x 



(2)\2 



,(4) 

'1111 



ox 1122 



24(4 3 2 ) 3 ) 2 



r(2) 



where 47 is the quadrupolar interaction and T the tem- 
perature in unergy units (ks — 1)- The quantities x^ 2 \ 
x^ and x^ are single particle thermal expectation val- 
ues [23[ , which are calculated by means of the single par- 
ticle potential. Numerical values of relevant parameters 
are given in Table 1. 

In case of a compressible lattice [41| the Hamiltonian 
H^ e L \ Eq. (l2"Tj) , gives the free energy term 



Fqqt/N — —2akp 2 [e x 



(26) 



where e xx etc. are the lattice strains. Since A < 0, this 
term favors e„ < 0, i.e. a volume contraction. The lattice 
Hamiltonian H leads to the clastic contribution 



Ftt/N = ^[c^« 



2 

XX 



L X X 1 



^yz 



(27) 



where V c = a 3 /4 is the volume per atom. 



In presence of an applied external pressure P, we con- 
sider the Gibbs free energy [4l| 



= FqQ + Fqqt + Ftt + PV(e 

xx T" & 



.'/.'/ 



(28) 



Minimizing G with respect to tu for i — x,y and z, we 
get an isostructural contraction of the cubic lattice 



txx = £ yy = e zz = - [8a 2 \A\p 2 + P]k l . 



(29) 



Here kl = (c^ + 2ci 2 )~ 1 is the compressibility. Isostruc- 
tural contraction of the lattice means that there is no 
symmetry breaking associated with the change of the cen- 
ter of mass positions of the Ce atoms at the 7 — > a tran- 
sition. We completely agree with conclusions from high- 
pressure and high-temperature x-ray diffraction experi- 
ments that the structure remains fee across the 7 — ^ tv 
transformation and retains crystallographic orientation 
during the transformation (42l |. Within the present the- 
ory a symmetry change occurs solely in the quadrupolar 
electronic charge densities whereby the Pa3 ordering of 
the latter is compatible with the fee structure of the lat- 
tice. 

Eliminating en from Eq. (|28|) and retaining only linear 
terms in P we rewrite G as 



G[p]/N = F /N + A'p 2 + Bp 3 + C'p 4 



where 



with 



A' = 



K (2) 



47 



(25) and where 



7 = 7 + a|A|«; L P, 



C = C -24a~ 1 A 2 n L , 



(30) 



(31) 



(32) 



(33) 



with C > 0. The occurrence of the cubic invariant in 
the free energy implies that the phase transition is of 
first order. The coexistence condition B 2 — AA'C = 
establishes a relation between the transition temperature 
T\ and the corresponding pressure Pi: 



Ti = x (2) (Ti)[47- 



B 2 {T X ) 
6C'(Tx) 



4a|A|K L P! 



(34) 



In first approximation Ti increases linearly with pressure, 
in agreement with the experimental (T, P) phase diagram 
of Ce jl| . Notice that the analogue of Eq. d34| also holds 
for Ceo fullerite [4l|, where experiments |43l, |44| show a 
linear pressure dependence of T% . 

With the parameters of Table Q] we obtain T\ = 89 
K. In accordance with the experimental situation [2TL 
[221 . [33l [HI we consider a fixed temperature T — T\ and 
a variable pressure. Then P< = P x + AP, AP < 0, 
corresponds to the disordered phase (in casu 7-Ce), while 



6 



AP (kbar) 

FIG. 2: Pressure dependence of the order parameter ampli- 
tude p, Pi taken as origin. 



P > = Pi + AP, AP > corresponds to the ordered 
phase (a-Ce). Minimization of G at T\ and P = P y 
with respect to p leads to the order parameter amplitude 



p(Ti,P > ) = 



8C'(Ti) 



Since B > 0, the order parameter p has to be negative 
and hence the — sign has to be chosen in front of the 
square root in Eq. (|35|) . The discontinuity of the order 
parameter on the transition line, Eq. (|34|) . is given by 



p(T 1 ,P 1 ) = - 



B(Ti) 
2C'{TiY 



(36) 



where we remind that T\ is an implicit function of Pi. 
With the parameters of Table Q] we obtain p\ = —0.0425. 
Here and in the following we assume that the quantities 
B(Tx), C'{Tx) and x (2) (Ti), all which are single expec- 
tation values, are not affected by moderate variations of 
pressure near Pi. In the ordered phase we obtain by 
means of Eqs. (133)) and (|34l) : 

p(Ti,P > )= 

-3B(T X ) - y/W(ft) 



192a|A|K L C"(Ti)AP 



8C'(Ti) 



(37) 



where AP > 0. Since p(T\, Pi) < 0, the order parameter 
amplitude increases in absolute value with the increment 
of pressure above Pi . In Fig. [5] we have calculated the 
order parameter amplitude by means of Eq. (|37[) . Here 
and below the cubic lattice constant a stands for a a in 
the ordered phase and a 7 in the disordered phase. 

In the following we will need the wave number depen- 
dent order parameter susceptibilities Xaa(k) m the disor- 
dered phase and in the ordered phase. The susceptibility 
is defined in terms of the fluctuations 8p( a \k) of the local 
order parameter from its average value: 



(38) 



Sp^(k) = p^(k)-VNpS m 



where i = x,y, z for a — 3,1,2, respectively. Hence, 



(39) 



We first consider the disordered phase. In case of a 
rigid lattice we obtain within molecular field theory 



1 



1 -T + x^vQe) 



(40) 



We recall that v(k) is the quadrupolar interaction matrix 
Eq. (IT0|) while 1 is the 3x3 unit matrix. The single par- 
ticle expectation value x^ = (p^ (n)p^ (ft)} has the 
same value for a = 1,2, and 3. Since we consider order 
parameter fluctuations near the phase transition, only 
wave vectors k near the star of k x are relevant. In ac- 
cordance with the condensation scheme (|lla[) . (|llbj) the 
matrix v(k) becomes diagonal at k = k x + 5k, where 
\Sk\ < \k*\. Expansion about k x gives 



-3P(Ti) - ^/9B 2 (Ti) ~32A'(Ti,P>)C'(Ti) 



v 22 (k)=-4 1+ ^-(5kl + 5kl) 



(41) 



(35 



imilar expressions with cyclic permutation of indices are 
obtained for vn(k) and ^(fc) near ky and fcf , respec- 
tively. In case of a non-rigid lattice we see from Eq. 
(f3"2"f that the applied pressure acts as addition to the 
quadrupolar interaction at k x . Replacing in Eq. (1411) 47 
by 47, we rewrite for the order parameter susceptibility 
in the disordered phase at temperature Ti and pressure 
P<: 



x p 2 p 2 (k) 



ie d (Ti,p<) 



a^[l + e d {Ti,P<){5kl + 5kl)\ 



Here the correlation length ^ is defined by 



^m,p<) 



x^(Ti) ia 2 



,1/2 



2(Ti-4 7 (P<) a; ( 2 )(Ti)) 



(42) 



(43) 



Expressions similar to Eq. (|4"2"|) are obtained for Xii(k) 
and X33(k) near k = k x and k x , respectively, and with 
(5kl + Ski) replaced by (<5fc 2 + 5k 2 z ) and (Ski + Sk^), 
respectively. 

Making use of Eqs. (f3~4"|) and (|3^|) . we obtain the cor- 
relation length 



3ja 2 C'(Ti) 



B 2 (Ti) + 24C'(Ti)a\A\K L AP] 



1/2 



(44) 



where AP < 0. With increasing pressure P < — > Pi in 
the disordered phase, the correlation length reaches its 
maximum value 



UTi,Pi) 



3 7 a 2 C"(Ti 



P 2 (Pi) 



1/2 



(45) 
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toward the onset of the phase transition. Likewise the 
susceptibility becomes maximum &tk = k* and (Ti,Pi): 



X22^x ) 



2 

a z 7 



(46) 



In the ordered phase the order parameter susceptibility 
is given by 



xZ(k) = j; {(p { Q\k) P %\-k)) - N P 2 % f 



(47) 



Applying mean-field theory [45[ we obtain for the 
quadrupolar model on a rigid lattice: 




x pp (fc) = X W U.T + x^v(k) + 1 • x^D(T)p 2 



(48) 



where 



D(T) 



T 



i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i 

-5 -4 -3 -2 -1 1 2 3 4 5 

AP (kbar) 

FIG. 3: Pressure dependence of the squared correlation 
length £ 2 (Ti,P) (in units A 2 ), Pi taken as origin. 



2(x( 2 )) 4 
5(^ } ) 



(2)^2 (4) o (4) 

X lm ^1122 



( X 123) 2 
,(2) 



where 



(49) 



a{Ti) = 3D(Ti) — C'{T\). 



(54) 



Numerical calculation show that D(T) > 0. In expres- 
sion (|47p p is the order parameter amplitude. Proceeding 
now as before we consider a nonrigid lattice and an ex- 
ternal pressure. Using Eqs. (|32|) and (|4Tj) we define the 
correlation length in the ordered phase at temperature 
Xi and pressure P > : 



£ Q (T 1( P>) = 



2(T 1 -4j(P>)xW(T 1 ) + xW(T 1 )D(T 1 )p2(T 1 ,P>)) 



Here the subscript Q stands for ordered quadrupolar 
phase. The corresponding order parameter susceptibil- 
ity for the component a — 2 is then given by 



Numerical evaluation of the quadrupolar model shows 
that a(Ti) > 0. Hence the correlation length decreases 
with increasing pressure in the ordered phase. Using 
Eq. flU) for £ d (Ti,P<) and Eq. ® for £ Q {T 1 ,P>) 
we have calculated the squared correlation length 
£c2(Q)(^l' P >< " < ^) as a function of pressure. The plot is 
shown in Fig. H From Eqs. 02J and (J3TJ) for X220) we 
see that Fig. [3] also describes the pressure dependence of 
the order parameter susceptibility at k = k£ ■ 

The expressions of the order parameter susceptibility 
(50) wm ^ e used in the next section where we study the tem- 
perature and pressure dependence of the elastic response 
near the first order phase transition. 



1/2 



ze (Ti,p' 



a 2 7 [l 



e Q (TuP>)(5kl 



■5kl)\ 



(51) 



Notice that the discontinuity of the order parameter at 
the first order phase transition leads to a drop of the 
correlation length and hence of the order parameter sus- 
ceptibility. With p(Tx,Pi) given by Eq. (j3"6"1) . we obtain 



(q(Ti,Pi) 



3 7 « 2 e'ero 



S 2 (T 1 )[1 + 3 J D(T 1 )/2C"(T 1 )] 



1/2 



IV. ELASTIC RESPONSE 

We will derive the static displacement - displacement 
response function matrix [46[ by using the well known 
Green's functions techniques [47j. Inversion of the re- 
sponse function in the long wavelength limit leads to the 
elastic constants. 

The Fourier transform of the retarded Green's function 
of two operators A and B is defined by 



(52) 

and since D/C > 0, £q(T x ,Pi) < ^ d (T u P{). The corre- 
lation length in the ordered phase is calculated by means 
of Eqs. ((50 ]) . ([37 ]) and ([51 ]) . The result reads 

Cq(Tx,P>) = {^C'iT,)) 1 ' 2 

[B 2 {T 1 ) + QC\T 1 )D{T 1 )p 2 {T ll P 1 ) 

+24aa(T 1 )|A|K L AP]~ 1/2 (53) 



((A;B)} 



dte lzt e(t)([A(t\B(0)}), (55) 



with frequency z = u) + ie, e — > + . The skew brack- 
ets (...) stand for a thermal average with the system's 
Hamiltonian H. In Heisenberg representation the time 
dependence reads A(t) = e mt Ae~ mt , H = H/h. We 
quote the equation of motion [47| : 



z((A;B)) z = ([A,B}) + (([A,n];B)) z 



(56) 
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with the identity 

(({A,H];B)) Z = -((A; {B,H})) 2 



(57) 



The dynamic displacement - displacement Green's func- 
tion is defined by 



D lJ (q,z) = ({in(q);uj(-q))) t 



(58) 



with lattice displacements and conjugate momenta speci- 
fied in Eqs. (I14a|i . (|14b|) . We recall that the Hamiltonian 
is given by Eq. (J]), where the parts have been speci- 
fied subsequently in Sect. [TTJ In the following we treat 
£l^(q) as a dynamic variable. Thereby we retain only 
the fluctuation contribution of the order parameter vari- 
able, writing 



(59) 



Applying twice the equation of motion (|56[) to Dij (q, : 



we use the commutation rules Eq. (|T5|) as well as the fact 
that the electronic variables p^ and QS % > commute with 
the lattice variables u and p. The result reads 



(z 6ik- M ik (q))D k j(q,z) 
h& ^ 



(60) 



Here the summation is understood over the repeated in- 
dex k. We recall that fl^(q) has been defined by Eq. 
(|20| . In a similar way we obtain by using Eq. (IBTl) 



(z 2 S tk - M ik (q))((n^(q);u k (-q))) z = 
lXa -. qj ((nM(q);rt j H-<f)))> 



(61) 



Taking the static limit z = in Eqs. and (I5T1) and 
combining the results we find 



A j (?> = 0) = Mr 1 ( 9 -) 
A 2 a 2 



—hS r 



Nm 



q r q, 



Mr 1 (5)«n«(g);nW(-5)»o 



Here the subscript stands for z = 0. 

We introduce the static susceptibilities [If 



-Dij(q,z = 0), 



X™(q) = -] i ((rt r) (q);rt s H-q)))o, 



(62) 

(63) 
(64) 



and rewrite Eq. (1B2"|) as 



Xl u (q) = M- 1 (q) 



AV 

Nm 



qrqsM-'^XrPiq) 



(65) 



We remind that in the long wavelength limit the static 
displacement susceptibility is related to the elastic con- 
stants [Ia| djki by 



(66) 



where p = 4m /a 3 is the mass density. Hence Eq. (|65|) 
allows us to calculate the elastic constants in presence of 
the coupling H^ e,L "> between electronic and lattice degrees 
of freedom. We recall that in absence of the coupling the 
bare elastic constants are given by Mi k (q) = p ~ qjqiC®j kl 
[see Eq. (|16|) ]. For cubic crystals we have in Voigt's no- 
tation the elastic constants C xxxx = en, C xyxy = C44, 
C xxyv = C12. With q = qi = g x (l,0,0) and i = j = x, 
Eq. (|65|) reduces in the long wavelength limit to 



Cll 



A 2 
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N 



X XX (9i 



0) 



(67) 



Here we have made use of Eqs. (J66f and (|16p . On 
the basis of the Hamiltonian, Eq. ((T|), the obtained ex- 
pressions for Xi^iq) an d en are rigorous results. Since 
X xx (qi = 0) > O7 the quadrupolar-elastic coupling, Eq. 
(IT^l) . leads to a reduction of the elastic constant en in 
comparison with the bare quantity c^. This reduction 
as a consequence of the ij( e>L ) coupling is responsible for 
the relative softening of the corresponding L[001] phonon 
branch in 7-Ce (32l. l35l| . In contradistinction to en, the 
shear elastic constant C44 is not affected by the coupling 
of the lattice to quadrupolar electronic density fluctua- 
tions. As we have shown in Sec. II, there is no coupling 
to shear strains e xy in H^ e ' L \ Eq. (fT9| . Hence taking 
A = in Eq. (|65| . we obtain for q = q\ and i = j = y: 



(68) 



(Xyy(qi)) 1 = M TO (gl), 



or equivalently by means of Eqs. (fl6|) and (|66| C44 = C44. 
Pressure experiments exhibit a discontinuity of the shear 
modulus [33j at the phase transition. However, this effect 
is solely due to the lattice contraction and a concomitant 
change of interatomic forces. Due to the absence of a 
direct coupling between lattice shears and orientational 
density fluctuations there are no precursor effects. 

We want to study the anomalous behavior of the elastic 
constant en near the quadrupolar ordering transition. 
From Eq. (|6"T)) it follows that the important quantity is 
the four-point (four factors p) susceptibility x X x (<fi = 0) • 
At high temperature and in the long wavelength regime 
such that the length scale of order parameter fluctuations 
is large in comparison with the interatomic spacing we 
use classical statistical mechanics. We then have 

X^(qi) = i<*fi (l) (ft)*fi (a) (-gi)>, (69) 



where 



SSlWfa) = rt x \qi) - (^(q!)). (70) 
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In Appendix we show that in the limit <fi — > 

! (<f 1 =0) = 2T^^'( X ^(fc)) 2 , (71) 



Xxx 



with a = 2, 3. Here ^(fc) is the order parameter sus- 
ceptibility for the component a. Cubic symmetry implies 
that 



£(*§£(*)) = E(x§5(*) 



(72) 



We calculate the right hand side of Eq. (|7ip by using 



E(^(^)) 2 



r 



(2tt) 3 



d 3 k{ X Z{k))\ 



(73) 



where = Na 3 /A is the volume of the crystal and N 
the number of atoms. So far the considerations of the 
present section are general and hold for the disordered 
as well as for the ordered phase. Since the expressions of 
the order parameter susceptibility in the disordered and 
in the ordered phase, Eqs. (|42|) and (1BT1) . respectively, 
exhibit the same wave vector dependence, we can treat 
both cases simultaneously. Since the integrand in Eq. 
(|73p does not depend on k y , the integration over dky in 
the interval ±27r/a yields 4ir/a. On the other hand the 
integrand vanishes for large values of |<5fc 2 + 8k^.\. Hence 
we integrate over a circle and extend the radius to oo. 
Taking into account Eq. (f72")l , we obtain 



no 

Xxx 



(9 = 0) 



2TA£ 2 

2 2~~ 5 



(74) 



where stands for ^ in the disordered phase and for 
£q in the ordered quadrupolar phase. Substitution of the 
result into Eq. (|67|) gives the longitudinal elastic constant 



Cll =c n [l + 3„(Ti,P)]- 



where we have defined 



^(ii,^J - -Q 



S 9 ' 



(75) 



(76) 



with = £d, Eq. (|4"4")l , in the disordered phase and £„ = 
£q, Eq. (j53")l , in the ordered phase. 

Near the phase transition in the disordered phase we 
have 

c 11 (T 1 ,F < )| d = C ; i [l + S d (T 1 ,P<)]- 1 . (77) 

The increase of £d(7i, P<) or equivalently of Ed(Ti,P < ) 
with AP < leads to a decrease of c\\ which reaches its 
minimum value 

cuCTx^OI, = c° n [1 + 3^(71, POP 1 (78) 

at P< = Pi. 



On the other hand at the onset of the phase transition 
one has 



Cll (T 1) p 1 )| Q = c? 1 [i + s Q (r 1 ,p 1 )] 1 



(79) 



The discontinuity of the order parameter and the con- 
comitant drop of the correlation length £q(Ti,Pi) < 
£d(Ti,Pi) results in a positive jump of en at the first 
order phase transition: 



cuCTi.POIq- e n (T 1 ,P 1 )\ d = 

ttoV [1 + SdCTi, Pi)][l + Sq(T x , J\)] ' 



(80) 



Here we have approximated a in the prefactor by (a a + 
a 7 )/2. 

In the ordred phase the decrease of the correlation 
length ^q(Ti,P > ) with increasing pressure leads to an 
increase of 



c n (T 1 ,P>)\ Q = c° n [1+H Q (T 1 ,P>)]" 



(81) 



The scenario described by Eqs. (l77|) - (|8Tj) is in full agree- 
ment with experiment. Indeed the pioneering exper- 
iments on elastic properties of Ce under pressure by 
Voronov et al. [34j as well as recent high resolution ul- 
trasonic measurements [HJ show that the propagation 
velocity of longitudinal ultrasonic waves decreases when 
the phase transition toward the a phase is approached 
with increasing pressure in the 7 phase. At the first order 
phase transition the longitudinal sound velocity exhibits 
a stepwise increase. With increasing pressure in the a 
phase the longitudinal sound velocity increases continu- 
ously. A corresponding behavior of the bulk modulus as 
a function of pressure was deduced from high-resolution 
neutron and synchrotron x-ray powder diffraction (2l| ex- 
periments. It was argued that the softening of the bulk 
modulus in the 7 phase with increasing pressure P — > Pi 
is a direct consequence of the softening of c\\. 

We close this section by observing that in Ceo fullerite 
there occurs a marked lowering of the bulk modulus B 
if at fixed T the fee— >sc phase transition is approached 
with increasing pressure (37j . The phenomenon is at- 
tributed to orientational reordering of the Cgo molecules. 
Given the similarities of the lattice related phenomena at 
the 7 — > a phase transition in Ce and the orientational 
phase transition in Cgo fullerite, we conclude that these 
transitions are isomorphic. From the mathematical point 
of view all essential aspects (Fm3m — > Pa3) are the 
same, while the constituents (electronic quadrupolar den- 
sities versus icosahedral molecules) are different. 



V. CONCLUSIONS 

The elastic properties of a model of interacting 4f- 
electron quadrupolar densities on a compressible fee lat- 
tice have been investigated as a function of temperature 
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I - H 




AP (kbar) 

FIG. 4: Pressure dependence of en (A = 5Ao), -Pi taken as 
origin. 



and pressure. On the basis of previous theoretical work 
[23l - [2a | this model, supported by recent nuclear spec- 
troscopy experiments [3l|, is taken as representative for 
the 7 — a "isostructural" phase transition in Ce. In par- 
ticular, we have studied by analytical theory the pres- 
sure dependent anomalies of the elastic constant en at 
the phase transition from the quadrupolar orientation- 
ally disoredered phase which we identify with the 7-phase 
(space group Fm3m) to the quadrupolar orientationally 
ordered phase which we identify with the a-phase (space 
group Pa3). As a result we find that en (equivalently 
the longitudinal sound velocity) decreases by approach- 
ing the phase transition from the disoredered phase with 
increasing pressure P < — > P\, at the first order phase 
transition en exhibits a positive jump, in the ordered 
phase en increases continuously with pressure P > > Pi. 
On the other hand we find that the elastic constant C44 
(equivalently the shear sound velocity) exhibits no pre- 
cursor effects near the transition. These theoretical re- 
sults are in full qualitative agreement with experiments 
[34,21,22,33] on elastic anomalies at the 7 — a "isostruc- 
tural" phase transition in solid Ce. 

We notice that the Hamiltonian H^ e,L \ Eq. (|17[) . which 
is quadratic in the electronic quadrupolar order param- 
eter variables and linear in the lattice displacements, 
accounts as well for the elastic anomalies as for the 
"isostructural" lattice contraction. The Hamiltonian is 
reminiscent from the compressible Ising model [49l |50T ] 
which is quadratic in the spin variables and linear in the 
lattice displacements. However, in the present case sym- 
metry properties are more subtle and account for the 
interplay between antiferroquadrupolar order in the elec- 
tronic densities and "isostructural" lattice contraction. 

Here a remark on symmetry reduction at the phase 
transition is in order. The quadrupolar interaction has 
its origin in the repulsive Coulomb interaction between 
4f-electrons on neighboring atoms. The ordering of the 
electronic quadrupoles on four sc sublattices reduces the 
repulsion and acts as an effective attraction. In reciprocal 



space the quadrupolar interaction matrix v(q), Eq. (|10p . 
has negative eigenvalues at the AT-point of the Brillouin 
zone of the fee lattice. 

Phonon dispersions of Ce measured by synchrotron ra- 
diation show pressure dependent anomalies related to the 
7 — a transition at the X-point of the Brillouin zone [32j . 
The explanation of these experiments is an outstanding 
challenge for theoretical work. 

A further problem is the quantitative improvement of 
the theory which should lead to increase the magnitude 
of the coupling coefficient A in Eq. (TT^j) for H^ e ' L \ At 
present we believe that the theory should be extended by 
including the intrasite coupling between 4f-electrons and 
conduction electrons (say 5d) and a possible influence of 
this coupling on the interatomic bonding. An indication 
of the relevance of such a mechanism is provided by recent 
work on the Ce dimer [l2j]. Likely such an extension of 
the theory would also contribute to elucidate the concept 
of a Mott transition 0, 0, HH versus a Kondo scenario 
at the 7 — a transition. 
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Appendix A 

Here we express (<Zi = 0) as a function of the order 
parameter susceptibility x^(fc). Starting from Eqs. (jol?)) 
and (|70[) we first notice that 



(fl (,) ®)) = (fi w («l = 0)) %,o 



Mi 



(Al) 



where 



(rt x \qi = 0)) = £'E<^ (a) (*)^ (a) ("*)>• (A2) 



We rewrite Eq. as 

X™(qi) = \ [<0<%-i)fi(*) (-«?!)) - (fi(*>(gi = 0)) 2 



(A3) 



We then approximate the four-point function 



<n«(ft)nM(-ft)> s £EEE 

k P a P 
(Sp^H-k + qi) 5p^ a \k) SpW\-f- qi) SpWtf)). 

(A4) 
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by the factorization scheme Substituting this result into Eq. (| A3|) and using Eq. ([39)1 

we obtain 

(fi<">(gi)n (B) (-«i)> * 

E E E'E' ^p [a) (-*) ^ (Q) (*)> ^ (/3) (0» 

fc P Q /3 

+ ( (-fc + &) 5pW (-p- £)) (5p(«) (fe) 8pW (p)) 

+(5pM {-% + q x ) &pW (p)} (8pW (k) SpW (-P- qi)) 

(A5) 



(9i) = 2r EE'^(fc-9i)x^W- (A7) 



or equivalently 

(QW(gl)QW(- g -l)) = (n<*>($i = 0)) 2 

+ 2 E E' (<^ (a) ^ (Q) (*-&)) 

x( ( 5 / 9 (a) (fc)(5 / 9 (a) (-fc)). (A6) Taking the limit £ = leads to the result of Eq. ([71]), 



[1] D.C. Koskenmaki and K.A. Gschneidner, Jr., Handbook 
on the Physics and Chemistry of Rare Earths, ed. K.A. 
Gschneidner, Jr., and L. Eyring (Amsterdam: North- 
Holland, 1978), p. 337. 

[2] P. W. Bridgman, Proc. Am. Acad. Arts Sci. 76, 55 
(1948). 

[3] M.R. MacPherson, G.E. Everett, D. Wohlleben, M.B. 

Maple, Phys. Rev. Lett. 26, 20 (1971). 
[4] L. D. Landau and E. M. Lifshitz.: Statistical Physics, 

(Pergamon, Bristol, 1995), Vol. 5. 
[5] A. W. Lawson, T.-Y. Tang, Phys. Rev. 76, 301 (1949). 
[6] B. Johansson, Phil. Mag. 30, 469 (1974). 
[7] B. Johansson, LA. Abrikosov, M. Alden, A.V. Ruban, 

and H.L. Skriver, Phys. Rev. Lett. 74, 2335 (1995). 
[8] A. P. Murani, S.J. Levett, and J.W. Taylor, Phys. Rev. 

Lett. 95, 256403 (2005). 
[9] A. Hjelm, J. Trygg, O. Eriksson, B. Johansson, and J. 
Wills, Phys. Rev. B 50, 4332 (1994). 
[10] X. Cao and M. Dolg, Mol. Phys. 101, 1967 (2003). 
[11] B.O. Roos, R. Lindh, P.-A. Malmqvist, V. Veryazov, and 

P.-O. Widmark, J. Phys. Chem. A 112, 11431 (2008). 
[12] A.V. Nikolaev, Phys. Rev. A 84, 012512 (2011). 
[13] J.W. Allen, R.M. Martin, Phys. Rev. Lett. 49, 1106 
(1982); J.W. Allen and L.Z. Liu, Phys. Rev. B 46, 5047 
(1992). 

[14] M. Lavagna, C. Lacroix, and M. Cyrot, Phys. Lett. A 
90, 210 (1982). 

[15] P. Fulde, Electron Correlations in Molecules and Solids 

(Springer Heidelberg, 1995). 
[16] L.Z. Liu, J.W. Allen, O. Gunnarsson, N.E. Christensen, 

and O.K. Andersen, Phys. Rev. B 45, 8934 (1991). 
[17] B. Johansson, A. V. Ruban, and I. A. Abrikosov, Phys. 

Rev. Lett. 102, 189601 (2009). 
[18] B. Amadon, S. Biermann, A. Georges, and F. Aryaseti- 

awan, Phys. Rev. Lett. 96, 066402 (2006). 
[19] J.-P. Rueff, J.-P. Itie, M. Taguchi, C.F. Hague, J.-M. 

Mariot, R. Delaunay, J.-P. Kappber, and N. Jaonen, 



Phys. Rev. Lett. 96, 237403 (2006). 
[20] G.Eliashberg, H.Capellmann, JETP Lett. 67, 125 (1998). 
[21] I-K. Jeong, TW. Darling, M. J. Graf, Th. Proffen, R. 

H. Heffher, Y. Lee, T. Vogt, J. D. Jorgensen, Phys. Rev. 

Lett. 92, 105702 (2004). 
[22] M.J. Lipp, D. Jackson, H. Cynn, C. Aracne, W. J. Evans, 

and A. K. McMahan, Phys. Rev. Lett. 101, 165703 

(2008). 

[23] A.V. Nikolaev and K.H. Michel, Eur. Phys. J. B 9, 619 

(1999) . 

[24] A.V. Nikolaev and K.H. Michel, Eur. Phys. J. B 17, 15 

(2000) . 

[25] A.V. Nikolaev and K.H. Michel, Phys. Rev. B 66, 054103 
(2002). 

[26] R. Sachidanandam and A.B. Harris, Phys. Rev. Lett. 67, 
1467 (1991); P.A. Heiney, J.E. Fischer, A.R. McGhie, 
W.J. Romanow, A.M. Denenstein, J. P. McCauly Jr., 
A.B. Smith HI, and D.E. Cox, ibid. 67, 1468 (1991). 

[27] W.I.F. David, R.M. Ibberson, J.C. Matthewman, K. 
Prassides, T.J.S. Dennis, J. P. Hare, H.W. Kroto, D.R.M. 
Walton, Nature, 353, 147 (1991). 

[28] W.I.F. David, R.M. Ibberson, T.J.S. Dennis, J.P. Hare, 
K. Prassides, Europhys. Lett. 18, 219 (1992). 

[29] P.A. Heiney, G.B.M. Vanghan, J.E. Fischer, N. Coustcl, 
D.E. Cox, J.R.D. Copley, D.A. Neumann, W.A. Kami- 
takahara, K.M. Creegan, D.M. Cox, J.P. McCauley, Jr., 
and A.B. Smith III, Phys. Rev. B 45, 4544 (1992). 

[30] F. Decremps, L. Belhadi, D.L. Farber, K.T. Moore, F. 
Occelli, M. Ganthier, A. Polian, D. Antonangeli, C. M. 
Aracne-Ruddle, and B. Amadon, Phys. Rev. Lett. 106, 
065701 (2011). 

[31] A.V. Tsvyashchenko, A.V. Nikolaev, A.I. Velichkov, 
A.V. Salamatin, L.N. Fomicheva, G.K. Ryasny, A. A. 
Sorokin, O.I. Kochetov, M. Budzynski, and K.H. Michel, 
Phys. Rev. B 82, 092102 (2010); A.V. Tsvyashchenko, 
A.V. Nikolaev, A.I. Velichkov, A.V. Salamatin, L.N. 
Fomicheva, G.K. Ryasny, A. A. Sorokin, O.I. Kochetov 



12 



and M. Budzynski, JETP 111, 627 (2010). 
[32] M. Krisch, D.L. Farber, R. Xu, D. Antonangeli, C. M. 

Aracne, A. Beraud, T.-C. Chiang, J. Zarestky, D. Y. 

Kim, E. I. Isaev, R. Ahuja, and B. Johansson, PNAS 

(USA) 108, 9343 (2011). 
[33] F. Decremps, D. Antonangeli, B. Amadon, and G. 

Schmerber, Phys. Rev. B 80, 132103 (2009). 
[34] F.F.Voronov, L.F.Vereshchagin and V.A.Goncharova, 

Doklad. Akad. Nauk SSSR, 135, 1104 (1960), Engl. 

transl. Soviet Phys. Doklady, 135, 1280 (1960). 
[35] C. Stassis, T. Gould, O.D. McMasters, and K.A. Gschnei- 

dner, Jr., and R.M. Nicklow, Phys. Rev. B 19, 5746 

(1979). 

[36] W. Schranz, A. Fuith, P. Dolinar, H. Warhanek, M. 
Haluska, and H. Kuzmany, Phys. Rev. Lett. 71, 1561 
(1993). 

[37] A. Lundin and B. Sundqvist, Phys. Rev. B 53, 8329 
(1996). 

[38] C.J. Bradley and A.P. Cracknell, The Mathematical The- 
ory of Symmetry in Solids, (Clarendon, Oxford, 1972). 

[39] P. Zielinski and K. Parlinski, J. Phys. C: Solid State 
Phys. 17, 3287 (1984). 

[40] K.H. Michel, J.R.D. Copley, D.A. Neumann, Phys. Rev. 
Lett. 68, 2929 (1992). 

[41] D. Lamoen, K.H. Michel, Phys. Rev. B 48, 807 (1993). 



[42] K.T. Moore, L. Belhadi, F. Decremps, D.L. Faber, J.A. 

Bradley, F. Occelli, M. Gauthier, A. Polian, and CM. 

Aracne-Ruddle, Acta Materialia 59, 6007 (2011). 
[43] G.A. Samara, J.E. Schirber, B. Morosin, L.V. Hansen, 

D. Loy, and A.P. Sylwester, Phys. Rev. Lett. 67, 3136 

(1991). 

[44] G. Kriza, J.-C. Ameline, D. Jerome, A. Dworkin, H. 
Swarc, C. Fabre, D. Schiitz, A. Rassat, P. BBernier, and 
A. Zahab, J. Phys. I (France) 1, 1361 (1991). 

[45] P.M. Chaikin and T.C. Lubensky, Principles of Con- 
densed Matter Physics, Cambridge University Press, 
Cambridge (1995). 

[46] W. Gotze and K.H. Michel, in Dynamical Properties of 
Solids, G.K. Horton and A. A. Maradudin Eds., North- 
Holland Publ. Comp. 1974, Ch. 9, p. 499. 

[47] D.N. Zubarev, Soviet Physics Uspekhi 3, 320 (1960). 

[48] The Princeton Companion to Mathematics, T. Gowers, 
Editor, Princeton University Press (2008), p. 27. 

[49] H. Wagner and J. Swift, Z. Phys. 239, 182 (1970). 

[50] D.J. Bergman and B.I. Halperin, Phys. Rev. B 13, 2145 
(1976). 

[51] K. Gofryk, E.K.H. Salje, D.J. Safarik, J.L. Smith, X. 
Ding, A.C. Lawson, P.S. Riseborough, B. Mihaila, and 
J.C. Lashley (submitted). 



